Influence of temperature correlations on phase dynamics and kinetics 
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The melting of ultrathin lubricant film is studied by friction between atomically flat surfaces. 
The fluctuations of lubricant temperature are taken into account denned by Ornstein-Uhlenbeck 
process. The phase diagrams and portraits are calculated for the cases of second-order and first- 
order transitions — the melting of amorphous and crystalline lubricants, ft is shown that in the 
first case the stick-slip friction domain appears, dividing the regions of dry and sliding friction. In 
the second case the three stick-slip friction domains arise characterized by transitions between dry, 
metastable and stable sliding friction. The increase of correlation time of the lubricant temperature 
fluctuations leads to increasing of the frictional surfaces temperature needed for realization of sliding 
friction. The singular point, meeting the mode of dry friction, has complex character of stability. 
The stationary states, corresponding to the stable and metastable sliding friction, are presented by 
the focus-type singular points. 

PACS numbers: 64.60.-i, 05.10.Gg, 62.20.Qp, 68.60.-p 
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I. INTRODUCTION 

The problems of the sliding friction continue to attract 
a considerable attention due to vast application in engi- 
neering of flat solid surfaces divided by thin film of lubri- 
cant It is found experimentally that during the 
process of friction the liquid film becomes progressively 
thinner, at first its physical properties change gradually 
(quantitatively), and then the changes acquire the sharp 
(qualitative) character. The qualitative changes consist 
of non-Newtonian mechanism of flow and replacement of 
the ordinary melting by glass transition, however, the 
lubricant film continues to behave itself like liquid. In 
tribology such behaviour is called the "mixed lubrica- 
tion" , which represents the intermediate mode of friction 
characterized by transition from volume properties of lu- 
bricating material to the boundary ones. 

The boundary mode of friction, described in the pro- 
posed work, is realized in the case of ultrathin lubricant 
films with thickness less than four diameters of molecules 
at smooth surfaces or asperities, high loads, and low shear 
rates. It is characterized by the following changes of 
static (equilibrium) and dynamic properties of lubricant 
— simple unstructured Newtonian liquid [3(: 

• non-fluidlikc (non-Newtonian) properties: transi- 
tion between liquid and solid phases, appearance 
of new liquid-crystalline states, epitaxially induced 
long-range ordering; 

• tribological properties: absence of flow until yield 
point or critical shear stress reached, solidlike be- 
haviour of lubricant characterized by defect diffu- 
sion and dislocation motion, shear melting, bound- 
ary lubrication. 
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Experiments with mica, silica, metal oxide, and sur- 
factant monolayer surfaces, between, which organic liq- 
uids and aqueous solutions were placed, have shown that 
there are transformations between the different types of 
dynamic phases during sliding They manifest them- 
selves in appearance of intermittent (stick-slip) friction 
HI IE IE 01 j which is characterized by periodic transitions 
between two or more dynamic states during the station- 
ary sliding and is the major reason for destruction and 
wear of frictional elements. Thus, molecularly thin lubri- 
cant films undergo more than one type of transition that 
results to existence of different types of stick-slip motion. 

In particular, the experimental investigation of Theo- 
logical properties of lubricant film were allowed to ex- 
plain the nature of above described anomalous features 
of boundary friction Besides, with this aim the the- 
oretical presentation of lubricating material was used as 
the viscoelastic matter characterized by heat conductiv- 
ity. So, the observed phenomenology of ultrathin lubri- 
cant film is successfully described in on the basis of 
viscoelastic matter approximation and Ginzburg-Landau 
equation, where order parameter defines the shear melt- 
ing and freezing. The phase diagram is calculated defin- 
ing the domains of sliding, stick-slip, and dry friction in 
the plane temperature — film thickness. 

Along the line 0, within the framework of Lorenz 
model for description of viscoelastic matter p| the transi- 
tion of ultrathin lubricant film from solidlike into liquid- 
like state is shown to occur due to both thermodynamic 
and shear melting. The cooperative description of these 
processes is carried out as a result of self-organization of 
the shear stress and stain, and the lubricant temperature. 
The additive noises of above quantities are introduced 
for building the phase diagrams, where fluctuations in- 
tensities and frictional surfaces temperature define the 
domains of sliding, stick-slip, and dry friction J3J. In ref- 
erence [lfj the conditions are found at which the stick-slip 
friction regime corresponds to the intermittency mode in- 
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herent in self-organized criticality phenomenon. 

However, in spite of the fact that, as a rule, temper- 
ature determines the lubricant state, the question re- 
mains opened on influence of its correlated fluctuations 
on the process of friction. Here, it is shown that inter- 
nal fluctuations of temperature, which are described by 
the Ornstein-Uhlenbeck process, result in complication 
of dynamic phase diagram. The shear stress distribution 
function, the phase diagrams and portraits are calculated 
for the cases of second-order and first-order transitions 
- the melting of amorphous and crystalline lubricant 
fSecs. iHl and IIIl)) . respectively. The different stick-slip 
friction modes are predicted to occur. The influence of 
correlation time value of the studied fluctuations is de- 
scribed on the phase diagram. The investigation of ki- 
netic modes of boundary friction is carried out using the 
phase plane method. The self-similar phase dynamics 
and kinetics of lubricant film are investigated. 



II. CONTINUOUS TRANSITION 

In the previous work [8j on the basis of rheological 
description of the viscoelastic medium the system of ki- 
netic equations has been obtained, which define the mu- 
tually coordinated evolution of the shear components of 
the stress a and the strain s, and the temperature T in ul- 
trathin lubricant film during friction between atomically 
flat mica surfaces. Let us write these equations using the 
measure units 

/ pc v rj T c \ 1/2 

e. = ^(^^(^Zi) 172 , Tc (1) 
Go \tt J \ Vo J 

for variables a, s, T, respectively, where p is the mass 
density, c v is the specific heat capacity, T c is the criti- 
cal temperature, 770 = rj{T = 2T C ) is the typical value of 
shear viscosity 77, tt = pl 2 c v /n is the time of heat con- 
ductivity, I is the scale of heat conductivity, k is the heat 
conductivity constant, r e is the relaxation time of matter 
strain, Go = i] Q /T e : 

T a a = -a + ge, (2) 
tJ = -e + (T - l)a, (3) 
t t T = (T e -T)-ae + cT 2 + \(t). (4) 

Here the stress relaxation time T a , the temperature T e of 
atomically flat mica friction surfaces, and the constant 
g = G/Go are introduced, where G is the lubricant shear 
modulus. Replacement of ejr a by de/dt reduces Eq. 
to the Maxwell-type equation for viscoelastic matter ap- 
proximation that is widely used in the theory of boundary 
friction The relaxation behaviour of viscoelastic lu- 
bricant during the process of friction is described also by 
Kelvin- Voigt equation © [1113 

It takes into account 



the dependence of the shear viscosity on the dimension- 
less temperature 77 = rj /(T — 1) 0. Equations J3J and 
(PJ) represent jointly the new rheological model. It is 
worth noting that rheological properties of lubricant film 
are investigated experimentally for construction of phase 
diagram Equation 10} represents the heat conduc- 
tivity expression, which describes the heat transfer from 
the friction surfaces to the layer of lubricant, the effect of 
the dissipative heating of a viscous liquid flowing under 
the action of the stress, and the reversible mechanic-and- 
caloric effect in linear approximation. Equations Q 
0) coincide with the synergetic Lorenz system formally 
[l3L Il4| , where the shear stress acts as the order parame- 
ter, the conjugate field is reduced to the shear strain, and 
the temperature is the control parameter. As is known 
this system can be used for description of the thermody- 
namic phase and the kinetic transitions. 

The purpose of this work is to study the phase dy- 
namics and kinetics of boundary friction at introduction 
into Eq. (0} the stochastic source X(t) representing the 
Ornstein-Uhlenbeck process: 

<A(t)) = 0, (A(i)A(t')) = - exp (Jl^l) , (5) 

where / is the temperature fluctuations intensity, t\ is 
the time of their correlation. 

Following ^] , let us define the physical sense of value 
/. The time correlation of dimensionless temperature of 
lubricant is determined by the average value of product 

tp(T) = (AT(t)AT(t + T)), (6) 

where AT(t) is the difference between the current and 
the average temperature. Within the framework of as- 
sumption on the quasistationarity of fluctuations AT the 
function ip(r) accepts the form: 

^(T)H(AT) 2 )exp(-C|T|). (7) 

Here 1 /Q fixes the relaxation time for setting of equilib- 
rium. Supposing that in formulas © and Q) the mo- 
ments of time are connected by equality t + r = t' , we 
get: 

<p(H-t) = (AT(t)AT(t')) = <(AT) 2 )exp(-C|<'-<|). (8) 

Then, substitution of average square of temperature fluc- 
tuations ((AT) 2 ) = T 2 /c v in expression © gives 

(AT(f)AT(f')) = — exp(-C|i' - t\). (9) 

Cy 

Comparing this formula with the second equality JSJ , we 
obtain: 

X(t) = AT(t), t x = - I = — (10) 

Consequently, the value of noise intensity / is defined by 
the temperature and the heat capacity c v of lubricant. 
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At the first glance it can appear, that in the phase dia- 
grams T e (I) presented further, the every value / meets 
the unique T e , and instead of regions, there the curve 
has to exist in pointed out coordinates. The motion 
along this curve describes the evolution of the system. 
However, it does not so, since parameter T e represents 
the thermostat temperature whose arbitrary change does 
not determine the temperature T uniquely. It is neces- 
sary only to suppose that as a result of system's self- 
organization the defined value T e can meet the manifold 
of values T and, in accordance with HlOfl . the intensities 
/ in the different moments of time. Besides, the variation 
of c v leads to the change of I in the course of time. Thus, 
it implies the existence of phase diagram. It is possible 
also to change the intensity I due to the arbitrary choice 
of parameter £ characterizing the concrete system. The 
last means that the systems exist, where the noise does 
not influence substantially on their time behaviour, and 
the ones, where fluctuations influence critically. 

However, here we do not restrict ourselves by such ap- 
proach and understand temperature fluctuations in more 
wide sense. It is related to that, as a rule, the thermal in- 
fluence is caused by different external stochastic sources. 
Besides, it is possible to interpret the noise with the help 
of fluctuations T, as well as of T e . This presentation of 
noise is often used for its modeling in real systems [l6[ . 

In reference || a melting of ultrathin lubricant film by 
friction between atomically flat mica surfaces has been 
represented as a result of action of spontaneously appear- 
ing shear stress leading to the plastic flow. This is caused 
by the heating of friction surfaces above the critical value 
T c o = 1 + g^ 1 . Thus, according to such approach the 
studied solidlike-liquidlike transition of lubricant film oc- 
curs due to both thermodynamic and shear melting. The 
initial reason for this self-organization process is the pos- 
itive feedback of T and a on e [see Eq. ©] conditioned by 
the temperature dependence of the shear viscosity lead- 
ing to its divergence. On the other hand, the negative 
feedback of a and e on T in Eq. plays an important 
role since it ensures the system stability. 

According to this approach the lubricant represents a 
strongly viscous liquid that can behave itself similar to 
the solid — has a high effective viscosity and still ex- 
hibits a yield stress 0, [llj . Its solidlike state corresponds 
to the shear stress a = because Eq. (J2Jl, describing 
the elastic properties at steady state & = 0, falls out 
of consideration. Equation containing the viscous 
stress, reduces to the Debye law describing the rapid re- 
laxation of the shear strain during the microscopic time 
t £ rs a/c r~j 10~ 12 s, where a ~ 1 nm is the lattice con- 
stant or the intermolecular distance and c ~ 10 3 m/s 
is the sound velocity. At the same time the heat con- 
ductivity equation takes on the form of simplest ex- 
pression for temperature relaxation that does not con- 
tain the terms representing the dissipative heating and 
the mechanic-and-caloric effect of a viscous liquid. 

Equation J3J describes the flow of lubricant with veloc- 
ity V = Ide/dt due to action of appearing viscous shear 



stress. Moreover, in accordance with Ref. 7] in the ab- 
sence of shear deformations the temperature mean-square 
displacement is defined by equality (u 2 ) = T/Ga. The 
average shear displacement is found from the relation- 
ship (u 2 ) — a 2 a 2 /G 2 . The total mean-square displace- 
ment represents the sum of these expressions provided 
that the thermal fluctuations and the stress are inde- 
pendent. Above implies that the transition of lubricant 
from solidlike to fluidlike state is induced both by heat- 
ing and under influence of stress generated by solid sur- 
faces at friction. This agrees with examination of solid 
state instability within the framework of shear and dy- 
namic disorder-driven melting representation in absence 
of thermal fluctuations. It is assumed that the film be- 
comes more liquidlike and the friction force decreases 
with the temperature growth due to decreasing activa- 
tion energy barrier to molecular hops. Besides, the fric- 
tion force decreases with increasing velocity at the con- 
tact V = Ide/dt because the latter leads to the growth 
of the shear stress a according to the Maxwell stress - 
strain e relationship: da/dt= —<r/T a +Gde/dt. 

This work is devoted to study of the stochastic source 
X(t) influence on the evolution of the stress <r(t). In ac- 
cordance with experimental data for the organic lubri- 
cant 0, 0] the stress relaxation time at normal pressure 
is equal to r a ~ 10~ 10 s. Since ultrathin lubricant film 
consists of less than four molecular layers, the temper- 
ature relaxes to the value T e during the time satisfying 
the inequality tt C t„. Therefore, we will suppose that 
conditions are fulfilled 

r ff «T e >r T , (11) 

at which lubricant temperature T follows the change of 
the shear components of stress a and strain s. Then, it is 
possible to select a small parameter and to put ttT ~ 
in Eq. (0J. As a result, we obtain the expression for 
temperature: 

T = T e -oE + o 2 + A(i). (12) 

Let us give for the system ©, ©, and (fT*S|) the more 
simple form, reducing it to the single equation for the 
shear stress a(t). For this purpose it is necessary to ex- 
press e and T via a. Differentiating with respect to time 
the equation for strain e that is obtained from J2J), we 
get equation for i. Substitution of these expressions for 
e, i and the equality l|12l) in Q gives evolution equation 
in canonical form of equation for nonlinear stochastic os- 
cillator of the van der Pole generator type: 

ma + 7 (a)a = f(a)+(l>(a)X{t), (13) 

where the coefficient of friction 7, the force /, the am- 
plitude of noise (f>, and the parameter m are defined by 
expressions 

7(d) EE - [T £ + T a (l + <T 2 )] , 

9 

f(a) ee a (T e - 1 - g- 1 ) - a 3 (g- 1 - l) , 

0(a) ee a, m=—. (14) 
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Let us find the distribution function of the stress a. With 
this aim, we will use the method of effective potential 01 j 
[Isj . As a result, the Fokker-Planck equation is obtained 
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It is expressed in terms of coefficients 

„-i 



^-M'-'*^(t + '-)}-< w > 

(17) 
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In the stationary case the solution of Eq. (|15|) leads to 
the distribution 



P(o-) = ,£ _1 exp{-£(cr)} 
which is fixed by effective potential 



E{a) = - J 



L>( 2 )(a 



da; 



(18) 



(19) 



and normalization constant 



(20) 



The stationary shear stress is found from the extremum 
condition of distribution (|TB)l 



£>W((x) 
£( 2 )(ct) 



= 0. 



(21) 



According to Fig. 1 the distribution 1|18|) has pronounced 
maximums whose positions are determined by the set of 
parameters r e , T\, g, I, and T e . At the small values of 
friction surfaces temperature T e a single maximum is re- 
alized at the point a = meeting the dry friction mode 
(curve 1). With T e growth the two maximums appear 
at points a — and a ^ (curve 2), the first of them 
corresponds to the dry friction, the second one — to the 
sliding. Here the stick-slip friction mode, characterized 
by transitions between the indicated stationary regimes, 
is realized. With the further growth of T e the zero max- 
imum of P{a) disappears and the maximum at a ^ 
remains only (curve 3), i.e., lubricant becomes liquidlike. 

Supposing in (|21() a = 0, we find the critical value of 
friction surfaces temperature 



led 
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2.9 



I 



(22) 



providing transition to the sliding mode. It is seen that 
T e0 grows at the increase of noise intensity / and cor- 
relation time t\. The values of relaxation times of the 



shear stress and strain influence by reverse manner. The 
domains of dry (DF), sliding (SF), and stick-slip (SS) 
modes of friction are realized in the phase diagram pre- 
sented in Fig. 2. The growth of correlation time t\ results 
in the increase of value T e , corresponding to the tricrit- 
ical point T at defined intensity I. Thus, the region of 
dry friction broadens, while sliding and stick-slip friction 
becomes more hardly realized. 

Apparently, the increase of the sheared surfaces tem- 
perature T e transforms lubricant to the sliding friction 
mode. It can be understood considering Eq. I|13[l that de- 
scribes the damping oscillations. Here the surfaces tem- 
perature is included only in expression for driving force 
/, which increases with growth of T e . As is known, the 
liquid state can correspond to the oscillation mode with 
large amplitude, but the solid state can not. With the 
increasing value of effective force in (|14(l the amplitude 
of oscillations increases, and longer oscillation process is 
realized to the moment of setting of the certain friction 
mode in the system. 

For studying the dynamics of change of friction modes 
it is enough to represent the distribution P(cr) by po- 
sition of its maximum a. This is achieved by the use 
of the path integrals formalism |l9j . within the frame- 
work of which the extreme values a — a(t) of the initial 
distribution function (|18f) evolve in accordance with the 
effective distribution 



n{<7, a} oc exp I — / A(<7, a, t)dt 



(23) 



Here the Onsager-Machlup function A, acting as the La- 
grangian of Euclidean field theory, is the subject for de- 
termination. 

Equation l|15fl can be transformed to the form: 



dP 
~dt 



d_ 

da 



D (l) 



da 



P 



da 2 



D (2) P) .(24) 



For finding of the A(a,a,t) dependence we write down 
the differential Langevin equation 



da = [ D 



dP^ 2 
da 



dt + V2DWdW(t) (25) 



corresponding to the Fokker-Planck equation Q24|l 19}- 
Here the stochastic differential dIF(t) represents the 
Winner process with the properties of white noise: 



(dW(t)) = 0, ((dW(t)) 2 ) = dt. 



(26) 



The feature of stochastic equations is that differential 
dW(t) can not be obtained by the simple division of 
Eq. (2>J by V2DW. With this aim, it is necessary to 
pass from the random process a(t) to the white noise x(t) 

1/2 

related with the initial Jacobian dx/d<r = (2Z?( 2 )) 
Then, substitution of Eq. into the Ito stochastic dif- 
ferential 



da; 



da; 1 d 2 x , ,_. 9 

d^ + 2d^ (d(7) ' 



(27) 
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taking into account i|26|) . leads to the expression 

^2 



dcr \ dcr 



dcr 



(28) 



Here the terms are neglected whose order exceeds 
(dW(t)). After reverse transition from the white noise 
x(t) to the initial process a(t) the equality 

dW(t) _ d 2DW + (DW)' 



dt V2D& 2V2D&> 



(29) 



is obtained, where the stroke stands for differentiation 
with respect to a. Plugging this expression into Gaussian 



and comparing with l|23|) . we arrive at Lagrangian 



A 



1 a 



~2 



4L>( 2 ) 



U 



with potential energy 



U = - 



16L»( 2 ) 



(30) 



(31) 



It is substantial that such form of potential energy U 
does not coincide with the effective potential . Con- 
sequently, for further consideration it is necessary to re- 
place the expression U <|31fl by E(a) l|19|l . In this case, 
Eq. (|30|l describes the system behaviour in accordance 
with distribution (|18|l obtained above. 

The system kinetics is defined by Euler-Lagrange 
equation 



d_8A 



dR 



(32) 



Within the white noise presentation the dissipative func- 
tion has the simplest form R = x 2 /2 and it is transformed 
to 



4L>(2) 



(33) 



with transition to the variable a — (2D^) 1 ^ 2 x. Substi- 
tuting in HT^i the equalities ll3l7|) . jJSt , t3lfl) . we arrive at 
differential equation 



. 2 p (2) y 

a 2Z?( 2 ) 



a — 9 r)( 2 ^ — 



(34) 



At the steady state supposing in 134|) a = 0, the ex- 
tremum condition (|21|l of distribution l|18l) is obtained. 
Note that the distribution maximum corresponds to the 
minimum of effective potential, and the distribution min- 
imum — to its maximum. 



Consider kinetics of the system using the phase plane 
method for analysis of Eq. I|34l) . The phase portraits are 
presented in Fig. 3 corresponding to the curves of Fig. 1. 

The region of dry friction (DF) (Fig. 3a) is charac- 
terized by the presence of the one singular point D, 
which corresponds to the maximum of probability P(cr) 
at a = 0. This point is non-standard and it requires in- 
terpretations. It is located at the origin of coordinates, 
and phase trajectories are curved around it, so that the 
system never comes to the value a — 0, i.e., this point 
is not stationary. Consider the system behaviour at the 
arbitrary initial condition. According to the phase tra- 
jectories the system evolves to the zero value of stress. 
Thus, if at initial conditions the growth rate of stress is 
positive it decreases to the zero at first (during this time 
the stress increases). Then, the stress value decreases 
asymptotically to the zero with the increase of its de- 
crease rate. Presumably, the described situation meets 
the mode of explosive amorphization, at which the sys- 
tem transforms very rapidly in amorphous solidlike state. 
The circumstance that zero stress is never achieved testi- 
fies to divergence of the probability P(u) at point a = 0. 
This is related with infinite growth of decrease rate of 
stress variation at going of the system to a = 0. Let us 
assume that the system has reached the point, at which 
decreasing rate becomes critical. In such case the in- 
crease of stress value is expected, which is accompanied 
by the change of rate sign and transition of the system 
in the positive phase plane domain. Further, again the 
decreasing of stress occurs, and the described situation 
repeats oneself. Alternatively, the sign of stress varia- 
tion rate does not change, and only its decreasing takes 
place. This moves the system on a neighboring phase tra- 
jectory, along which it approaches to the zero stress more 
rapidly. As a result, in the course of time the oscillation 
mode of dry friction is set in the vicinity of point a = 
at the arbitrary initial conditions. Thus, the oscillations 
amplitude is small, and lubricant has solidlike structure. 

The phase portrait of the system, characterizing the 
region of stick-slip friction (SS), is shown in Fig. 3b. The 
three singular points appear here: D, the saddle N , and 
the stable focus F. As well as above, the point D is re- 
alized at the origin of coordinates and it corresponds to 
the dry friction mode. The saddle N meets the mini- 
mum of P(cr) and it is unstable stationary point. It is 
worth noting that at the initial value of the shear stress 
on the right-hand side from point N and a = 0, the slid- 
ing mode of friction is set in the system during time. If 
the initial value of stress appears on the left-hand side 
from saddle N, the dry friction is set in similar case. 
Thus, the point N separates two maximums of distri- 
bution function P{cr). The focus F corresponds to the 
non-zero maximum of stress distribution function, i.e., 
it describes the liquidlike state of lubricant. The damp- 
ing oscillations, corresponding to this point, mean that 
lubricant becomes more liquidlike and more solidlike pe- 
riodically. However, the stable sliding friction is set al- 
ways. Presumably, these oscillations are conditioned by 
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the presence of noise. 

The phase portrait, meeting the sliding friction (SF) 
(Fig. 3c) characterized by single non-zero maximum of 
distribution function P(a), has one singular point — the 
stable focus F. The stability of sliding friction is con- 
firmed also by large overexpansion of phase trajectories 
near F along axes of ordinates and abscissas. However, 
it is apparent that at the initial large value a the system 
does not reach the point F, and approaches asymptoti- 
cally to the zero stress value. This circumstance implies 
that conditions can be realized, at which the system will 
be near the regime of dry friction. As described above, 
at reaching the value of critical variation rate of stress its 
sign changes and becomes positive. It is seen from the 
phase portrait that in such case the system also passes 
to the mode of stable sliding friction. 



III. INFLUENCE OF DEFORM ATIONAL 
DEFECT OF MODULUS 

Actually, the shear modulus introduced (in terms of 
the relaxation time t ct ) in Eq. iffil. depends on the stress 
value. This leads to the transition of the elastic defor- 
mation mode to the plastic one. It takes place at char- 
acteristic value of the stress a p , which does not exceed 
the value a s (in other case the plastic mode is not mani- 
fested). For consideration of deformational defect of the 
modulus we will use T a (a) dependence proposed in [8J, 
instead of r CT . As a result, Eq. J5J takes the form: 



£>M and D^: 



-all 



1 



1 + a/a 



+ 5ee, 



(35) 



where the relaxation time for the plastic mode t p = rj a /0 
is introduced {r] a = r a G is the effective viscosity, is 
the hardening factor), # = 0/G<lis the parameter 
describing the ratio of tilts of the plastic and the Hookean 
sections of the deformation curve, ge = G 2 /0Go and a = 
a p /a s are the constants. Then, within the framework of 
approximation pif) the system l|35|l . ©, and @, as well 
as above, is reduced to equation (cf. (|13fl ): 



ma + j(a)a = f(a) + 0(<r)A(i), 



(36) 



where the coefficient of friction 7, the force /, the am- 
plitude of noise <p, and the parameter m are defined by 
expressions 



7(0-) 



- 1 



(1 + a/af 



+ r p (l + a 2 



T.-l-i 



+ a/a 



9 V 1 + a/a 



+ a/a 



1+a/a 



(ff)=<r,m=SZ£. (37) 
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According to the effective potential method 01 1 \IM we 
obtain the Fokker-Planck equation (|15|l with coefficients 



- < a 

7 



g \ 1 + a /a 
_1 + a/a 



21a 



1 

1 - - . 

9 V 1 + a I a 

a ( tJI - a ~ l 
7 



Iar\ 



d& = r 7 -i 

7 



g \ (1 + a/a) 3 a 
2r A l . 



,(38) 
(39) 



Using Eqs. l|2T]l. and (|2HJl, the equality is found 

with the same meaning as (|22H , which gives the boundary 
of existence of distribution l|18fl maximum at zero value 
of stress corresponding to the solidlike state of lubricant 
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The P(a) dependencies are shown in Fig. 4 for the 
different modes of friction. Curves 1-5 correspond to the 
domains of dry (DF), stick-slip (SS), stick-slip and sliding 
(SS+SF), metastable and stable sliding (MSF+SF), and 
sliding friction (SF), respectively. The phase diagram 
and portraits are presented in Figs. 5 and 6 corresponding 
to the curves of Fig. 4. 

The most complex form of P(a) function is inherent 
in SS+SF region (curve 3 in Fig. 4). Here the solid- 
like, metastable and stable liquidlike lubricant states co- 
exist meeting the maximums of P(a). It means the real- 
ization possibility of intermittent (stick-slip) friction, at 
which the periodic transitions occur between the dynamic 
modes corresponding to these states. In the MSF+SF 
domain of stick-slip motion the metastable and stable 
sliding can periodically change each other (curve 4 in 
Fig. 4). It is characteristic that transition from SS+SF 
to MSF+SF region is accompanied by disappearance of 
dry friction in the system. The domain of dry friction 
(DF) broadens, and the region of sliding friction (SF) 
decreases with growth of correlation time t\ of noise. 

The phase portrait of dry friction region (DF) is similar 
to that is inherent in continuous transformation (Fig. 3a) . 
It implies that DF domain is equivalent at the taking into 
account of the modulus defect and without it. 

The phase portrait describing the region of stick-slip 
friction (SS) is similar to the characteristic one at contin- 
uous transformation (Fig. 3b). The basic their difference 
is that here the trajectories around focus are considerably 
more elongated along both coordinates axes. It means 
the greater stability of sliding friction. 

The most complex region (SS+SF) is represented by 
the phase portrait shown in Fig. 6a. The five singular 
points are realized here: D, the saddles N, N', the sta- 
ble focuses F, F' . As well as above, saddles correspond 
to the minimums of P{a) dependence. Point D meets 
the solidlike state of lubricant. Stable focus F determines 
the first non-zero maximum of probability. It is apparent 
that the oscillations are weakly pronounced around this 
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point. In this mode lubricant represents the very viscous 
liquid, because in such type of fluid at presence of noise 
the oscillations are damped strongly. Actually, the point 
F corresponds to the small values of stress, and with its 
decreasing the lubricant becomes more viscous, and it is 
transformed into the solidlike state at a = 0. Thus, using 
phase portraits it is possible to give explanation to that 
the liquidlike state of lubricant, but not the solidlike one, 
corresponds to the large values of shear stress. Focus F' 
meets the second non-zero maximum of P(<r) function, 
and there are oscillations with large amplitude around it. 
This implies the fluidlike state of lubricant and, accord- 
ingly, sliding. The last point is in large distance along 
abscissa axis from all others ones. This mode of fric- 
tion is most probable only, since points D and F, corre- 
sponding to the dry and metastable sliding friction, have 
large stability and probability of realization also. From 
here the conclusion follows that lubricant can undergo 
periodic transitions (stick-slip) between the modes corre- 
sponding to the points D, F, and F' . Since these modes 
are stable and separated by the pronounced minimums of 
distribution function P(a) (by saddles), the transitions 
between them is necessary to expect after large intervals 
of time. 

The phase portrait of MSF+SF domain is represented 
in Fig. 6b. There are three singular points — the stable 
focuses F, F' , and the saddle N. The latter is similar 
to the described above saddles and it meets the mini- 
mum of probability dependence on the stress. Point F 
corresponds to the first maximum of distribution, which 
describes metastable sliding mode (MSF), and F' - 
the second maximum, which defines stable sliding (SF). 
There are only insignificant oscillations around the focus 
F, however, lubricant in this mode is less viscous liquid 
than in the vicinity of the point F' in Fig. 6a. At the 
origin of coordinates the singular point is absent, and 
the dry friction is not realized. Focus F' is similar to 
described in Fig. 6a, however, it's "attraction" domain 
is more stretched along both axes that means the larger 
fluidity of lubricant and stability of this mode. Therefore 
in comparison with the previous case, here the arising of 
sliding friction is more probable. 

SF region is represented by the phase portrait, which 
is similar to the described for the continuous transforma- 
tion (Fig. 3c). Here the one stable focus F is realized, 
characterized by oscillations in it's vicinity, representing 
stable sliding friction. The basic difference is that in this 
case oscillations take place with large amplitude that im- 
plies strong fluidity of lubricant and pronounced stability 
of such mode. However, as well as in all above considered 
situations, in the course of time in lubricant the station- 
ary shear stress is set corresponding to the maximum of 
the initial distribution P{cr). 

In the basic equations © - J3J the shear stress a 
stands in the first power. However, in general case it's 
exponent a may be fractional (0 < a < I), but not inte- 



ger: 

T a a = -<r a +ge, (41) 

r £ i = -e+(T-l)a a , (42) 

t t T = {T e -T)-a a e + a 2a + X(t). (43) 

Taking into account the additive noises of the shear stress 
and strain, and the temperature of lubricant film it has 
been shown [ljj that such system describes the self- 
similar mode for which the characteristic scale of shear 
stress is absent pcj . This regime is determined by the 
homogeneous distribution function 

P(y)=y- 2a V(a), y = aa s . (44) 

In particular, the value 2a = 1.5 corresponds to the self- 
organized criticality mode, at which, unlike the phase 
transition, the process of self-organization does not re- 
quire the external influence (T e — 0) and occurs sponta- 
neously mim. 

The self-similar behaviour of lubricant film is studied 
taking into account fluctuations of its temperature de- 
fined by Ornstein-Uhlenbeck process. It is shown that 
the fluctuations of lubricant temperature result in dis- 
appearance of sliding friction region at presence of dry 
and stick-slip friction domains in both cases of second- 
order and first-order transitions. In the second case the 
stick-slip motion arises characterized by three stationary 
values of shear stresses at which dry, metastable and sta- 
ble sliding friction are realized. The increase of correla- 
tion time of lubricant temperature fluctuations leads to 
increasing of frictional surfaces temperature needed for 
realization of stick-slip friction. 

The study of Eqs. jjjl - shows that the phase 
portraits for similar domains reproduce above considered 
qualitatively. However, there are substantial differences. 
The fractional Lorenz system at a ^ 1 and I results 
in presence of the singular point D in phase portraits, 
which corresponds to the solidlike state of lubricant. Be- 
sides, the variation of a leads to the complication of P(a) 
dependence, and as a result, to more complex form of 
phase portraits. Within the determined friction mode 
at decreasing of a the increase of abscissas of the stable 
focuses is observed. Consequently, the weakening of frac- 
tional feedbacks in the Lorenz-type models results in the 
increase of lubricant fluidity and the reducing of friction. 
However, in the systems described by fractional exponent 
a the dry friction is realized always. Thus, it is impossi- 
ble to assert that such systems more preferable to friction 
decrease than linear systems. 

IV. CONCLUSION 

The above consideration shows that increase of tem- 
perature of frictional surfaces T e , at presence of colored 
noise of lubricant temperature, can be accompanied by 
self-organization of elastic and thermal fields leading to 



8 



the mode of sliding friction. Indeed, the correlation de- 
gree of lubricant temperature change plays the substan- 
tial role. If the correlation time t\ increases, the growth 
of friction surfaces temperature is necessary for transition 
from dry to sliding friction mode at the fixed intensity 
/ of temperature fluctuations. In the case of continu- 
ous transformation at small intensity I this transition 
takes place avoiding the region of intermittent friction, 
i.e., has the form of second-order transition — the melt- 
ing of amorphous lubricant. In the reverse case of large I 
the first-order transition is realized corresponding to the 
melting of crystalline lubricant. 

At setting of the sliding friction mode in the system the 
damping oscillations arise and the shear stress relaxes to 
the stationary value fixed by the probability distribution 
P{u). The amplitude of these oscillations increases with 
growth of stationary values of shear stress. It means that 
large shear stress a corresponds to the liquidlike structure 
of lubricant. The solidlike state of lubricant is described 
by the singular point D at the origin of coordinates that 
has complex character of stability and corresponds to the 
divergence of probability P{cr). The oscillations near this 
point are absent. 

For description of first-order transition the shear mod- 
ulus defect is taken into account. It is shown that the 
change of temperature fluctuations intensity / and fric- 
tional surfaces temperature T e can transform the system 
from the dry friction mode to the sliding. The latter 



arises at two values of shear stress. Accordingly, the three 
singular points, that define the stationary values of stress, 
appear in phase portraits — non-standard point D at 
zero stress and two stable focuses at non-zero ones. The 
intermittent (stick-slip) mode of friction can be realized 
as a result of transitions between solidlike, metastable 
and stable liquidlike lubricant states, which are described 
by zero and non-zero singular points. 

Taking into consideration the nonlinear relaxation of 
the shear stress and fractional feedbacks in the Lorenz 
system it has been shown that the region of sliding fric- 
tion is realized in the phase diagram only in absence of 
temperature fluctuations. In this case the singular point 
D is always present in phase portraits, which corresponds 
to the solidlike state of lubricant and dry friction. Be- 
sides, the change of fractional exponent can complicate 
the phase diagram, increasing the number of friction do- 
mains, in accordance with experimental data Q. 
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Figure captions 

to the paper by A. V. Khomenko 
"Influence of temperature correlations on phase 
dynamics and kinetics of ultrathin lubricant film" 

FIG. 1. The distribution function of the shear stress for 
the second-order transition at g = 0.2, t ct = t £ = 
0.1, t a = 0.2, and 7 = 5. The curves 1, 2, 3 cor- 
respond to the temperatures T e — 5, 16, 20, respec- 
tively. 

FIG. 2. The phase diagram corresponding to the pa- 
rameters of Fig. 1 with the domains of dry (DF), 
sliding (SF), and stick-slip (SS) friction (T is the 
tricritical point). 

FIG. 3. The phase portraits corresponding to the pa- 
rameters of Fig. 1: (a) DF mode corresponds to 
the curve 1 in Fig. 1; (b) SS — curve 2 in Fig. 1; 



(c) SF — curve 3 in Fig. 1. 

FIG. 4. The distribution function of shear stress for 
the first-order transition at t p — r e = t\ = 
0.1, 0- 1 = 7, a = 0.3, g e = 0.4,7 = 4.5. 
The curves 1-5 correspond to the temperatures 
T e =16, 21, 23.25, 24, 26, respectively. 

FIG. 5. The phase diagram corresponding to the pa- 
rameters of Fig. 4 with the domains of dry (DF), 
sliding (SF), and stick-slip (SS, MSF+SF, SS+SF) 
friction modes. 

FIG. 6. The phase portraits corresponding to the pa- 
rameters of Fig. 4: (a) SS+SF mode corresponds 
to the curve 3 in Fig. 4; (b) MSF+SF — curve 4 in 
Fig. 4. 
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